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Linear  spectral  coherence  me  assures  have  been  used  in  the  neuro¬ 
sciences  to  test  hypotheses  which  address  the  question  of  whether 
multiple  E0G  recording  sites  are  independent  or  whether  they  are 
activated  by  cannon  sources  of  neurophysiological  activity,  '’’his 
measure  is  appropriate  when  regional  neural  sources  interact  and  thereby 
electrically  activate  several  recording  sites  through  linear  trans¬ 
mission  pathways  which  may  or  may  not  be  different  in  their  linear 
transformation  properties.  However,  if  the  transmission  media  are 
nonlinear,  then  interactive  dependency  is  not  necessarily  revealed  by 
a  linear  coherence  test.  Therefore,  if  carman  sources  are  activating 


EEG  recording  sites  through  nonlinear  media,  evaluating  the  resulting 
relationships  among  the  signals  recorded  from  these  sites  requires  a 
test  which  reveals  the  presence  of  such  nonlinear  relationships.  In 
neurophysiological  applications,  a  polycoherence  cross-spectral  measure 
provides  such  a  test  far  nonlinear  dependency  (similar  to  the  linear 
coherence  test)  among  E3G  recording  sites.  The  data  requirements  and 
statistical  properties  of  these  linear  and  nonlinear  measures  are 
described  and  results  of  a  linear  coherence  analysis  are  presented  in 
the  context  of  an  EEu  pilot  study  of  learning-disabled  children. 


INTRODUCTION 


o 

Analysis  of  the  relationships  among  multiple  channel  EBG  signals  ' 
represents  one  of  the  challenging  problems  in  the  evaluation  of  brain 
electrical  activity.  One  frequently  used  analysis  procedure  for  the  pair-wise 
coirparison  of  multiple  channel  data  is  provided  by  the  cross-correlation 
function,  which  measures  the  product  moment  correlation  between  two  signals  as 
a  function  of  lag  (i.e. ,  the  time  displacement  between  the  signals).  However 
this  procedure  has  the  disadvantage  of  being  insensitive  to  the  presence  of 
highly  correlated  low  amplitude  frequency  components  of  activity  which  may  be 
imbedded  in  uncorrelated  high  amplitude  activity.  Therefore  the  normalized 
Fourier  transform  of  the  cross-correlation  function(a)  is  commonly  employed  to 
provide  a  means  of  deriving  the  degree  of  linear  dependence  or  coherence 
between  two  signals  over  distinct  frequency  bands.  This  allows  linearly 
dependent  activity  contained  within  particular  frequency  bands  to  be  clearly 
observed  in  the  presence  of  intense  uncorrelated  or  linearly  independent 
activity  in  other  frequency  bands. 

The  question  of  dependence  between  two  signals  can  be  broadened  to 
include  nonlinear  dependence.  However,  to  perform  an  analogous  test  for 
nonlinear  dependence  requires  the  use  of  polycoherence  analysis  which  involves 
the  conput&tion  of  higher  dimensional  Fourier  transforms.  This  paper  will 
describe  the  analytical  methodology  for  testing  of  both  linear  and  nonlinear 
dependence.  Illustrative  pilot  data  from  a  study  using  EBG  linear  spectral 
coherence  to  test  hypotheses  involving  shared  brain  activity  ir.  learning 
disabled  children  as  conpared  to  normal  controls  will  also  be  described. 

WThe  normalized  Fourier  transform  of  the  cross-correlation  function  is 
called  the  "spectral  coherence  function"  or  simply  the  "coherence  function". 
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The  issue  of  neurophysiological  interest  examined  in  this  paper  can  he 
stated  as  follows:  can  EEG  activity  recorded  from  multiple  scalp  locations 
give  information  about  the  level  of  shared  or  cannon  electrical  activity 
associated  with  the  brain  regions  proximate  to  these  recording  sites?  A  basic 
method  for  addressing  this  question  is  to  test  whether  a  transformation  exists 
which  will  map  signal  S\(t)  into  signal  S2(t)  where  S\(t)  and  S2(t)  represent 

EEG  signals  recorded  from  two  scalp  locations. 

If  the  mapping  can  be  dare  by  a  linear  transformation  then  the  signals  in 
question  are  said  to  be  linearly  dependent.  The  degree  of  linear  dependence 
(as  a  function  of  frequency)  between  the  ERG  signals  associated  with  these 
brain  regions  is  determined  by  applying  spectral  coherence  analysis. 

Similarly,  if  the  mapping  of  S\(t)  into  S2(t)  requires  a  nonlinear 
transformation  then  poly spectral  (in  particular  cross-bicoherence f  analysis 
can  be  used  to  address  the  question  of  coupling  or,  stated  another  way,  the 
degree  of  nonlinear  dependence  as  a  function  of  frequency. 

The  degree  of  dependence  among  signals  can  range  from  zero  for 
independent  signals  to  unity  for  totally  dependent  signals.  In  the  context 
of  EEG  analysis  this  degree  of  dependence  or  connectivity  is  taken  as  a 
measure  of  shared  electrophysiological  activity  among  brain  regions. 

LINEAR  DEPENDENCE 

The  numerical  measure  of  linear  dependence  is  derived  from  the  coherence 
function  whose  magnitude  ranges  from  zero  to  unity.  The  analytical  and 
statistical  methods  involved  in  applying  the  coherence  function  to  EEG  data 
are  outlined  in  the  following  paragraphs. 

Hie  coherence  function  is  defined  in  terms  of  the  normalized  cross- 
spectrum  of  the  two  time  series.  The  cross-spectrum  can  be  expres  xi  as 
the  Fourier  Transform  of  the  cross-correlation  function  as  follows: 
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Let  S\(t)i  S2(t)  denote  the  different  time  series 
and  let  $i2(t)  denote  the  cross-correlation  function  of  Si (t),  S2(t), 
where  x  is  the  variable  time  shift  between  S i  and  S2,  then  is 

defined  by  the  integral  equation 


$12^; 


Sx(t)  s2(t  +  \)dt 


(1) 


The  cross-spectrum  of  S\(t)  and  S2(t) ,  denoted  as  P\2(f)  is  defined  by 


(2),  using  the  exponential  form  of  the  Fourier  transformation, 


-  /  <tn2(i) 

•V.OO 


(where  /  is  frequency  in  Hertz  and  x  is  time  shift  in  seconds). 

By  substituting  equation  (1)  into  equation  (2)  and  appropriately  factoring 
the  resulting  double  integral  it  can  be  shown  that 


Pi  2  (f)  =  Si(f)  s$  (f)  (3) 

where:  Si(f)  is  the  Fourier  transform  of  S,  ( t)  and  S2(f)  is  the 
conjugate  Fourier  transform  of  S2(t).  Since  S\(f)  and  S*(f)  are  complex 
numbers  in  general,  therefore  P\2(f)  is  also  complex.  It  simplifies 
subsequent  analysis  to  represent  them  in  polar  form,  as  follows: 


Sx(f)  =  i>iezdl 

(4) 

S*2(f)  =  r2e"102 

(5) 

Substituting  (4)  and  (5)  into  (3)  gives  the  cross-spectral  density  in  polar 
form: 

Pi2(f)  =  rir2e^^1~Q2^  =  v\v2e^  (6) 

(where  0  =  0i  -  62) 

While  not  explicitly  shown,  note  that  the  r’ s  and  0's  are  functions  of 
f.  Now  to  arrive  at  the  spectral  coherence  function  we  normalize  equation 
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(6)  by  dividing  by  nr2  and  then  average  over  a  number  of  frequencies  and/or 
over  a  number  of  time  epochs  of  the  two  time  series  being  analyzed. 

Thus  if  the  average  is  taken  over  2N+1  discrete  frequencies  we  may  write 
the  spectral  coherence  as 


C<fn>  = 


0'  2N+1 


,Mfk> 


where  fQ  is  the  center  of  the  spectral  window  and  uniform  weighting  is  used 
in  the  window.  If  the  average  is  taken  over  M  time  epochs  at  a  specific 
frequency,  fv,  (ensemble  averaging)  we  may  write  the  spectral  coherence  as 


C(fk)  =~jr1L  eiQm(fk)  (8) 

m=l 

where  is  the  phase  angle  between  S\  and  S2  at  frequency  during 

the  mth  time  epoch.  If  a  cento inat ion  of  frequency  and  ensemble  averaging 
is  used  then  spectral  coherence  may  be  written 


C(fo]  *  M(2N+1) 


Again,  these  expressions  assume  the  use  of  uniform  weighting  in  the 
spectral  window. 

It  should  be  noted  that  far  linearly  independent  signals  the  phase 
difference  0(=Qi-92J  is  a  randan  function  over  both  frequency  and  ensemble 
and  therefore  the  expected  value  of  coherence  averaged  over  frequency  and/or 
time  (ensemble  averaging)  is  zero. 

Unis,  when  coherence  deviates  significantly  from  zero,  one  may  conclude 
that  the  two  signal  processes  in  question  are  related  through  a  linear  trans¬ 
formation  over  the  spectral  region  where  such  significant  deviation  occurs. 

The  statistical  measure  of  significance  of  the  deviation  of  a  coherence 
es+  '  .iate  fran  its  expected  value  depends  on  the  number  of  independent  samples 
of  coherence  which  are  used  in  arriving  at  the  coherence  estimate,  and  upon 


the  probability  distribution  function  of  these  coherence  estimates  as 
described  below. 


A  coherence  estimate  is  the  average  ccnputed  from  samples  of  the 
normalized  cross-spectrum.  These  sanples  can  be  represented  as  points  on 
the  unit  circle  in  the  complex  plane,  as  illustrated  in  Figure  1.  If 
these  points  are  uniformly  distributed  over  the  unit  circle,  it  is  clear 
that  the  expected  value  of  coherence  lies  at  the  origin  (x  =  y  =  0) .  The 
normalized  samples  of  cross-spectrum  may  be  written: 


_  e^k  =  cosQ.  +  i  sinQ-, 


(10) 


iyk  —  &  ~  —  UL/&  ~r  u  o  w* 

A  random  distribution  of  values  on  the  circumference  of  the  unit  circle 
(as  shown  in  Figure  1)  is  equivalent  to  the  condition  that  the  relative 
phase  values  between  Si  and  S2  be  uniformly  randomly  distributed  over  an 
ensemble  of  time  epochs  and/or  over  a  nurrber  of  frequencies  in  a  spectral 
window  centered  at  frequency  /  .  Under  this  condition  the  complex  value  of 
coherence  may  be  written: 


N  N 

Coh(f0>  -  -§•  Z  **  *  i  -y  Vk  (11) 

Since  x k  and  y k  lie  on  the  unit  circle 
x k  =  cosQ^ 
yk  =  sinQk 

the  PDF  (Probability  Distribution  Hinction)  of  both  x  and  y  (given  that  0 
is  a  uniformly  distributed  random  variable)  have  the  same  form  and  are  given 
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The  corresponding  means,  ( x ,  y)  and  variances  (a2,  a2)  are  given  h 

x  y 


x  =  y  =  0 


where  N  =  number  of  samples . 


(13) 

(14) 


Fran  (13)  and  (14)  the  resulting  mean  and  variance  of  the  coherence 
magnitude  (\aoh(f0)  |  5  v)  are: 

r  -  0  (15) 


Thus,  to  test  the  hypothesis  that  S\(t)  and  S2(t)  are  linearly 
independent  time  series,  one  examines  the  probability  that  the  empirical! y- 
obtained  mean,  computed  from  N  samples,  deviates  from  the  expected  value 
(zero  in  this  case). 

For  example  when  the  coherence  magnitude,  r,  is  obtained  by  a  combination 
of  ensemble  averaging  over  20  epochs  and  frequency  averaging  over  5 
spectral  components  then 


N  =  20  x  5  =  100 


and  the  resulting  standard  deviation  of  the  estimate  is 
1  1 

a  =  —  =  — - —  =  Q.10 

r  sr  mr 

Thus  an  empirically  obtained  value  of  coherence  magnitude  which  exceeds 
0.2  would  be  more  than  two  standard  deviations  from  the  expected  value  for 
independent  signals.  Therefore,  it  would  be  statistically  reasonable  to 
conclude  that  the  signals  in  question  are  not  independent,  but  rather  that 
they  are  linearly  dependent  or  coherent  over  the  spectral  region  where  the 
coherence  magnitude  exceeds  0.2. 

NONLINEAR  DEPENDENCE 

The  test  for  a  nonlinear  relationship  between  S\(t)  and  S2(t)  is 
similar  in  principle  to  the  linear  test,  but  it  makes  significantly  greater 
demands  on  ccnputer  processing. 
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By  analogy  with  the  linear  coherence  test,  the  nonlinear  test  involves 
a  milt i-dimentional  Fourier  transform  of  a  multiple  lag  product  similar  to 
the  cross-correlation  function,  whose  normalized  Fourier  transform  is  the? 
coherence  function.  Specifically  the  double  lag  product  is  given  by: 


«  -y 

/u\ 

The  cross-b ispect rum  '  ' 
of  <j)  (x iji2^  as  follows: 


Sift)  Sift  +  Tji  S2(1.  +  r2)  dt  (17) 

is  defined  as  the  two  dimensional  Fourier  transform 


Bff ■  f  )  =  f*°  f**  ^(tuX2)  e-j2*fmU  e-02*fni*d^dr. 

J—OO  J—CO 


(IB) 


The  cross-bicoherence  is  obtained  by  normalizing  the  cross-bispectrum  as 
follows: 

Cross-Bicoh  (f  3  f  )  =  Average  over  frequency  and/or  ensemble  of  the 
following  triple  product  of  F.T.’s 

Sjlfj  sllfn>  sUf^  +  4; 

\s,(fal\]sl(fn)\\st(fm  +  fj\ 


where  4  denotes  canplex  conjugate  and  the  vertical  bars  denote  absolute  value. 
It  is  instructive  to  write  this  triple  product  in  polar  form,  as  follows: 
Denote 


S'(V  •Tme 


% 


S*(f  +  f  )  =  r  e~^P 
1  Jm  Jn  p 


(19) 

(20) 

(21) 


where  the  r’s  and  <J>Ts  are  functions  of  frequency  and  \s^(f)\  =  ry  Using 
the  polar  form  in  equations  (19),  (20)  and  (21),  cross -bicoherence  is  given 
by 


Cross- Bicoh  (f^,  fn)  =  average  over  frequency  and/or  ensemble  of 


(by 


If  Sift )  =  Sift),  then  B  (fm. 


f  )  is  called  the  bispectrum. 


ft 


where  ^  =  <f>m  +  <pn  -  <pp  compares  the  phase  at  the  sum  frequency  ( fm  +  J'y) 

in  channel  2  with  the  sum  of  the  phases  at  frequencies  /  and  in  channel  1 . 

Thus  for  every  pair  of  frequencies  selected  we  obtain  a  value  for  f., 
which  defines  a  point  on  the  unit  circle  in  the  complex  plane  (just  as  in 
the  case  of  linear  coherence).  Therefore  if  ( 2N  +  1)  pairs  of  frequencies 
are  selected  over  an  ensarble  of  M  time  epochs,  the  bicoherence  estimate  is 
obtained  by  averaging  the  (P.N  +  1)M  samples  as  follows  (assuming  uniform 
weighting  of  the  samples) : 

1  M(  2N+1)  . 

Cross-Sicoh  (fm,  fn)  =  M-(2N  --  i}  £  e  (22) 

k= 7 

Some  detailed  expansion  is  required  to  implement  the  above  simplified 


expression .  First,  the  bicoherence  estimate  for  the  frequency  pair  ,  /  ) 
is  an  average  obtained  from  samples  centered  around  (/  ,  f  )  as  follows: 


Bicoh  (f  f)  = 


m*  J n  2N+1 


v  5'l'4  *  w  s*lfm  *  4  *  MW  (23) 

k~!l  |S,  ("4  *  Up  |  Is,  (fn  *  Uf!  I  |s2C4  +  4  +  sup  | 


Using  a  subscript  J  to  designate  the  time  epoch  over  which  the  above 
estimate  is  obtained  and  Af  to  designate  the  frequency  increments,  then  when 
ensemble  averaging  as  well  as  frequency  averaging  is  used  we  may  wife 


.  2  fr  si,/ D.V=up 

Btaoh  (f.  f)  ~  TrrrzjTTT  2-j  - — - (2d) 

*  k~N  \SUj(fm+kAf)  \  |  SUj(fn+Hf)  I  \S2,j(ft+fn+2kAf)  | 


v  'm 


The  argument  for  the  statistics  of  this  estimator  is  similar  to  the  previous 
discussion  of  linear  coherence.  Thus  for  two  independent  random  processes 
the  standard  deviation  (a)  of  bicoherence  magnitude  from  the  expected  value 


of  zero  is  o  = 


mW+T) 


Thus,  if  one  averages  5  pairs  of  lines  (N=2)  and  over  20  epoches  (M=20) 
then  a  =  0.1,  similar  to  the  example  for  linear  coherence. 
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OCHERENCE  ANALYSIS  OF  PILOT  DATA 

Linear  coherence  analysis  was  applied  to  pilot  data  fran  the  EEGs 
of  white  male  learning  disabled  children  and  normal  controls.  The 
objective  of  this  analysis  was  to  investigate  whether  differences  in 
EEG  multiple  channel  linear  dependence  provides  a  basis  for  separating 
these  populations. 

The  coherence  analysis  results  shewn  in  Figure  2  represent  TUT 

data  recorded  from  a  learning  disabled  child  and  an  age  and  sex  mi  ed 

normal  control.  The  recording  paradigm  was  designed  to  control  tf 

cognitive  state  of  the  subjects  so  that  data  could  be  obtained  dur.  t, 

information  processing.  Although  several  tasks  were  employed,  the 

typical  results  in  Figure  2  were  derived  fran  EEG  data  recorded  while 

the  subject  indicated  whether  sequentially  presented  pictures  and  words 

matched.  The  subject  was  instructed  to  delay  the  response  until  a  tone 

was  presented  approximately  three  seconds  after  the  second  stimulus  and 

then  press  a  button  if  the  stimuli  matched.  Hie  coherence  analysis  is 

based  on  the  1.2  seconds  of  EEG  in  the  middle  of  this  delay  interval. 

The  most  consistent  finding  which  separated  the  12  learning  disabled 

from  the  13  normal  children  in  the  pilot  study  was  coherence  magnitude 

at  frequencies  above  20  Hertz.  It  can  be  noted  in  Figure  2  that  for 

frequencies  above  20  Hertz  the  learning  disabled  child  exhibits  values 
(c) 

of  coherence  '  much  less  than  0.1  while  the  normal  control  exhibits 
values  of  coherence  greater  than  0.1  over  most  of  this  frequency  range. 

The  comparisons  in  Figure  2  are  for  coherence  between  left  and  right 
central  leads  in  the  test  and  control  subjects,  and  similarly  for  coherence 
between  left  and  right  occipital  leads. 

777^ - 

'As  in  the  discussion  of  linear  coherence  the  coherence  estimates  here 

are  based  on  averaging  over  20  epochs  and  5  spectral  lines.  Ebr 

linearly  independent  signals  the  standard  deviation  from  zero  of  these 

estimates  Is  0.1. 
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These  findings  suggest  that  impaired  sharing  of  high  frequency  EEF 
activity  may  be  implicated  in  the  neural  processes  associated  with 
learning  disability.  More  generally  these  spectral  procedures  for 
obtaining  dependency  measures  among  multiple  ERG  recording  sites  offer 
a  quantitative  basis  for  assessing  regional  electrophysiological 
interactions  during  mental  processing  or  in  response  to  medication  or 
other  clinical  intervention. 

OCNGLUDING  COMMENTS 

In  evaluating  the  EEGs  of  subjects  it  is  important  to  not  that 

major  differences  in  EEG  waveshape  across  subjects  may  be  of  no  eonse- 
3  4 

quence  ’  as  far  as  revealing  useful  diagnostic  or  clinical  differences. 
The  essential  property  of  multiple  channel  coherence  analysis  is  that  it 
measures  the  degree  of  camrunication  or  shared  EEG  activity  among  the 
channels  recorded  from  a  particular  subject  rather  than  waveshape  per  se. 
It  is  this  coherence  measure  of  shared  activity  which  constitutes  the 
criterion  for  comparing  or  categorizing  subjects.  Thus  even  though  the 
EEG  waveshapes  across  subjects  may  be  quite  different  it  is  possible 
that  the  multiple  channel  shared  EEG  activity  of  such  subjects  may  be 
quite  similar.  Ultimately  it  is  hoped  that  the  miltiple  channel  measures 
described  here  can  provide  important  information  ^o  aid  in  diagnosis  and 
In  the  evaluation  of  treatment  of  conditions  such  as  learning  disability 
as  well  as  provide  a  measure  of  mental  processing  or  neurophysiological 
state  which  is  independent  of  idiosyncratic  differences  in  EEG  waveshape. 
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